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ABSTRACT 

This paper presents numerical solutions of jet-induced mixing in a par- 
tially full cryogenic tank. An axisymmetric laminar jet is discharged from the 
central part of the tank bottom toward the liquid-vapor interface. Liquid is 
withdrawn at the same volume flow rate from the outer part of the tank. The 
jet is at a temperature lower than the interface, which is maintained at a cer- 
tain saturation temperature. The interface is assumed to be flat and shear- 
free and the condensation-induced velocity is assumed to be negligibly small 
compared with radial interface velocity. Finite-difference method is used to 
solve the nondimensional form of steady state continuity, momentum, and energy 
equations. Calculations are conducted for jet Reynolds numbers ranging from 
150 to 600 and Prandtl numbers ranging from 0.85 to 2.65. The effects of 
above stated parameters on the condensation Nusselt and Stanton numbers which 
characterize the steady-state interface condensation process are investigated. 
Detailed analysis to gain a better understanding of the fundamentals of fluid 
mixing and interface condensation is performed. 


NOMENCLATURE 


Aj surface area of central jet 

Aout surface area of outflow at the tank bottom 

A s surface area of the liquid-vapor interface 

B tank to jet diameter ratio, D/d 

Cl.C 2 .C 3 .C 4 constants 

specific heat at constant pressure 
tank diameter 
jet diameter 

gravitational acceleration 
condensation heat-transfer coefficient 
latent heat of condensation 


D 

d 

9 

h c 

h fg 


Ja b 

J Sj 

k 

m c 

m,n 

Nu c 

P 

Pg 

P* 

0 

Re j 

r 

r* 

St c 

T 

T j 

T S 

T* 

u 

u c 
u j 

u out 

u* 

u ; 

V 

v s 

V* 


bulk Jacob number, Cp(T s - Tout )//h fg 
jet Jacob number, C p (T s - Tj)/hfg 
thermal conductivity 
condensation mass flux 

constants of exponent in the correlation equations 

condensation Nusselt number, h c D/k 

pressure 

equilibrium hydrostatic pressure 

dimensionless pressure, (p - Pg)/pu? 

jet volume flow rate 

jet Reynolds number, puj d/p 

radial coordinate measured from the centerline 

dimensionless radial coordinate, r/D 

condensation Stanton number, h c /pUjCp 

temperature 

jet temperature 

interface temperature 

dimensionless temperature, (T - Tj)/(T S - Tj> 
axial velocity 

condensation-induced velocity 
jet velocity 
outflow velocity 

dimensionless axial velocity, u/uj 
dimensionless condensation-induced velocity, u c /uj 
radial velocity 

radial velocity at the interface 
dimensionless radial velocity, v/uj 
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v* dimensionless radial velocity at interface 

x axial coordinate measured from the tank bottom 

x* dimensionless axial coordinate, x/D 

Xp potential core length 

x s 1 iquid fi 1 1 ing height 

Greek symbols: 

ij dynamic viscosity 

p liquid density 

— average value over the interface 


Subscripts : 


j 

eval uated 

at 

out 

evaluated 

at 

s 

evaluated 

at 

c 

evaluated 

at 

ct 

eval uated 

at 

b 

eval uated 

at 


jet inlet 
outflow location 
liquid-vapor interface 
condensation condition 
tank centerline 
bulk liquid 


INTRODUCTION 

The thermal environment in space can result in radiation heat leak into 
the cryogen storage tanks and increase the tank pressure. Of the various tech- 
nologies (ref. 1) being developed for tank pressure control, axial jet-induced 
mixing (ref. 2) is considered as an efficient method for a short-term storage 
system. The mixing of the tank contents causes a forced convection and removes 
heat quickly from the interface into the bulk liquid. The induced interface 
condensation then results in the reduction of tank pressure. It is obvious 
that the vapor condensation on the subcooled liquid at the interface plays a 
key role in reducing the tank pressure. 

Several studies of steady-state fluid mixing and condensation rate have 
been conducted. Thomas (ref. 3) measured the condensation rate of steam on 
water surfaces mixed by a submerged jet. The liquid height to tank diameter 
ratio was varied from 0.32 to 1.3. It was found that the condensation rate was 
roughly proportional to the jet Reynolds number. Dominick (ref. 4) conducted 
experiments to investigate the effects of jet injection angle and jet flow rate 
on the condensation rate in a Freon 113 tank. The ratio of liquid height to 
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tank diameter was about one. The tests showed that heat transfer at the inter- 
face was enhanced by the mixing process as jet injection angle became more nor- 
mal to the interface. The average heat transfer coefficient at the interface 
was found to be increasing with jet Reynolds number to a power of 0.73. Sonin, 
Shimko, and Chun (ref. 5) also measured the steady-state condensation rate of 
steam in a water tank with liquid height to tank diameter ratio greater than 
three. The condensation process was dominated by the liquid-side turbulence 
near the interface. About 0.3 of tank diameter below the interface, the turbu- 
lence was found isotropic and essentially unaffected by the proximity of the 
surface. A correlation between the condensation rate and characteristic inter- 
face turbulent velocity was also developed in reference 5. Hasan and Lin 
(ref. 6) used a finite-difference method to numerically solve time-averaged 
conservation equations along with K-e turbulence model for the prediction of 
turbulent velocity. Since the information about the turbulence at the inter- 
face was not available, zero gradients of K and e were assumed for the 
boundary condition applied at the interface. The numerical prediction was in 
good agreement with Sonin's data except for the region close to the interface 
where diffusion process was dominant. The failure of numerical solution for 
the near-interface region, as stated in reference 6, was probably due to the 
inappropriate turbulence model or the inappropriate turbulence boundary condi- 
tions applied at the interface. 

For laminar flows, the fluid properties and boundary conditions are well 
defined. This motivates the present study to investigate more fundamentals of 
fluid mixing and the associated interface condensation process by a submerged 
laminar jet. For a typical cryogenic storage tank, the maximum liquid height 
to tank diameter ratio is about one and the mean-flow velocity is generally 
much larger than the turbulence at the interface. It is hoped that the study 
of laminar jet-induced mixing can provide some useful information for studying 
the turbulent jet-induced mixing in a typical cryogenic tank. In space, the 
effect of buoyancy force may be insignificant due to the low-gravity environ- 
ment. Therefore, the fluid mixing in low-g may be equivalent to the mixing 
of constant-density liquid in normal gravity provided that the effect of free 
surface configuration is excluded. The present study is conducted for a 
normal-g environment with the neglecting of buoyancy force due to constant 
density. Numerical solutions are obtained by using a finite-difference method 
to solve nondimensional continuity, momentum and energy equations. Calcula- 
tions encompass a wide range of the Prandtl numbers and jet Reynolds numbers. 
Detailed analyses are performed and simplified equations are generated to 
describe the effects of the above parameters on the condensation Nusselt and 
Stanton numbers which characterize the steady-state condensation process at 
the liquid-vapor interface. 


THE PHYSICAL PROBLEM CONSIDERED 

The physical system and the coordinates used to analyze the problem are 
shown in figure 1. A circular cylindrical tank of diameter D contains liquid 
with a filling height x s . An axisymmetric laminar jet with a velocity Uj 
is discharged from the central part of the tank bottom toward the liquid-vapor 
interface. The liquid with same volume flow rate, 0, is withdrawn from the 
outer part of the tank bottom such that the liquid fill level is kept constant. 
The outflow area, A ou t, is much larger than inflow jet area, Aj , (about 357:1) 
such that the outflow velocity is very small compared with the jet velocity. 
Both the central jet and outflow velocities are assumed uniform. The ratio of 
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jet diameter, d, to tank diameter, D, is fixed at 1:20. The ratio of liquid 
height to tank diameter, x s /D, is set to be one. 

The liquid-vapor interface is assumed to be flat (wave free) and shear 
free and at a constant saturation temperature, T s . All the solid walls are 
insulated. The central jet is kept at a constant temperature, Tj,, lower than 
the interface temperature. All the associated thermodynamics and transport 
properties are assumed constant since only small temperature variation is con- 
sidered in the whole flow field. The outflow region is assumed to have zero 
temperature gradient. 

The energy balance at the interface yields the following relation between 
condensation mass flux and heat transfer 


Ws 


■I 


m c h fg 


dA s - 




dA. 


( 1 ) 


where m c and m c are the average and local, respectively, condensation mass 
fluxes at the interface, hfg is the latent heat of condensation, and A s is 
the interface surface area. Experimentally, the steady-state condensation rate 
can be obtained based on the measurement of the temperature difference between 
outflow and inflow jet (ref. 5). The relation of energy balance for the whole 
system yields 


m c h fg A s = 


PQ C o (T out - V 


(1 


Ja b ) 


( 2 ) 


where Q is the jet volume flow rate, Ja^ is the bulk Jacob number defined as 
C p (T s - T ou t)/hfg, and T ou ^ is the outflow temperature which is essentially 
uniform over the outflow area. Equations (1) and (2) can be used as a checkup 
of the numerical solutions. For commonly used cryogens such as hydrogen and 
nitrogen, Ja^ is usually much less than one. For calculation convenience, 

Jat, is assumed zero in the present study. This assumption will be justified 
by the obtained solution later on. 


A local condensation heat-transfer coefficient, h c , and a local condensa- 
tion Stanton number, St c , which describe the interfacial heat and mass trans- 
ports can be defined as 


h 


c 


m cS 

<T S V 


(3) 


It is noted that, instead of the bulk liquid temperature and the interface 
velocity, the jet temperature and the jet velocity, respectively, are used in 
equations (3) and (4). There are two reasons for these unusual definitions. 
First, the jet conditions are generally easy to control and measure so that the 
definitions in equations (3) and (4) may be more useful. Second, the differ- 
ence between bulk temperature and interface temperature is very small in the 
present study so that h c will be very sensitive to the numerical error if 
(T s - Tb> is used in equation (3). For convenience, the average values of h c 
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and St c at the interface are generally used in the analysis and are expressed 
as 


r-^fL 

c <T s - V 


St. = 


c " pU 3 C P 


(5) 


( 6 ) 


MATHEMATICAL FORMULATION 

The jet-induced mixing considered in the present study is steady-state and 
incompressible with gravity acting in the vertical -x direction. The govern- 
ing equations for this axisymmetric laminar flow problem are: 

!_<„> „ JL (rv) , o 


fx Clj2) ♦ r37 <urv) 


1 8p 

P ax + 


g + 


y 9 2 u 
p 9x 2 



9_ (UV ) + _i_ (rv 2 ) = _ 1 9fi _ + k_9_v + IL JL ( r dv) 

dx (UV ' + r9r rv p 9r p f 2 + p a>( 2 + p r9r V 9 rj 


h < uT> + 4? 


k_ 9 h k 9 ( 9T\ 

p 9x 2 + pC p r9r V dr) 


(vrT) = r + r 

pC n 3v 2 P C 


The boundary conditions are applied for the system sketched in figure 1 
At the centerline, the symmetric conditions are used: 


v = o — = — = o 

u ’ 9r 3r 


At the solid walls, the nonslip and adiabatic conditions are employed: 


u = v = 0, 


8T 9T 
9x 9r 


0 


The free surface is assumed to be flat (wave free) and shear free. The 
condensation-induced velocity at the interface, u c , is equal to 


u 


c 



(7) 


and is assumed negligible compared with surface velocity. This assumption 
will also be justified by the obtained numerical solutions. The neglecting of 
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condensation-induced velocity u c leads to the uncoupling between fluid dynam- 
ics and thermal calculations. Thus, the free surface boundary conditions for 
velocity are. 


u = 0, g.O 


The temperature of the interface is kept constant at T s , 
the velocity and temperature are assumed to be uniform: 


For the jet inlet, 


u = Uj, 


v = 0, 


T = T- 


Liquid is withdrawn from the outer part of the tank bottom with the same volume 
flow rate as that of the injected liquid jet, so that the uniform withdrawing 
velocity can be obtained by Q/A ou t . The conditions for radial velocity and 
temperature applied at this liquid-withdrawing plane are given by 

Nondimensional ized variables can be formed as 


x* -i 


r* - 1 
r D’ 


u* = 

u . 
3 


v* = 

u . ’ 


p* = 


p - p. 


pu 


J 


T S T j 


where, for convenience, the gravity term has been subtracted from the 
x-momentum equation by using the static equilibrium equation: 


% 

9x = - pg 


The dimensionless forms of the steady state governing equations are 



9u* 9r*v* 



8x* + r*9r* 
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9u* 2 9u*r*v* 

+ 1 

9 2 U 

9x* + r*9r* " 

9x* + B Re. 

vj 
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r*9r 


8u* 

9r* 


9u*v* 

9r*v* 2 

9p* 1 v* 

1 

9 2 v* 9 

9x* 

+ r*9r* 

" 9r* ' B Re^ f *2 B Re^ 

[ 9 x* 2 + r * 9r * 


9u*T* 

9r*v*T* 1 

fg 2 t* 

3 / * 9T*' 


9x* + 
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J 
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The relevant parameters in the governing equations are the jet Reynolds number 
(Rej), the Prandtl number (Pr) and tank to jet diameter ratio (B): 



pU i d 
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C qP 

Pr - 
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The average condensation Nusselt number at the interface, Nu c , is defined by 


From equations (5) and (6), the 
number can then be expressed as 


- h c D 

Nu c = -V 


average condensation Nusselt number and Stanton 


Nu = 8 
c 


(IS) 


dr' 


s 

Nu" 

St - — 

ir c B Re, Pr 


The corresponding boundary conditions then become: 
At the centerline (symmetric conditions) 

3u* 3T* 


v* = 0, 


3r” 


3r 


* = 0 


( 8 ) 

(9) 


At the solid walls (nonslip and adiabatic conditions) 


At the interface 


At the jet inlet 


u* = v* = 0, 


3T* 3T* n 

3x* = 3r* = u 




v* = 0, 


T* = 1 
T* = 0 


At liquid-withdrawn plane 


u* = 


u j A out 


V* = 0, 



The effects of jet Reynolds number ( Re j ) on the dynamics of fluid mixing 
in the tank and the steady-state condensation rate at the liquid-vapor inter- 
face are investigated. Calculations will be performed with the jet Reynolds 
number ranging from 150 to 600 in which the flow is expected to be laminar 
according to the experiments of McNaughton and Sinclair (ref. 7). The Prandtl 
number (Pr) will be varied from 0.85 to 2.65 to represent various cryogens 
such as hydrogen and nitrogen so that the performance of fluid mixing can be 
evaluated. 
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NUMERICAL METHOD OF SOLUTIONS 


The above dimensionless forms of elliptic partial differential equations 
are numerically solved by a finite-difference method. The finite-difference 
equations are derived by integrating the differential equations over an elemen- 
tary control volume surrounding a grid node appropriate for each dependent 
variable (ref. 8 ). A staggered grid system is used such that the scalar prop- 
erties, p and T, are stored midway between the u and v velocity grid 
nodes. The bounded skew hybrid differencing (BSHD) is incorporated for the 
convective terms (ref. 8 ) and the integrated source terms are linearized. 
Pressures are obtained from a predictor-corrector procedure of the Pressure 
Implicit Split Operator (PISO) method (ref. 9) which yields the pressure 
change needed to procure velocity changes to satisfy mass continuity. The 
governing finite-difference equations are solved iteratively by the ADI method 
with under relaxation until the solutions are converged. 

Calculations are performed with a nonuniform grid distribution with con- 
centration of the grid nodes in the centerline, near-wall, and near-interface 
regions where the gradients of flow properties are expected to be large. The 
nonuniform grid distribution in axial direction is generated by using an expo- 
nential function of Roberts' transformation (ref. 10) with stretching parameter 
equal to 1 . 02 . In the radial direction, the scheme (ref. 11 ) with a constant 
ratio between two adjacent grid spacing is used. This ratio is set to be about 
1.22 in the present study. To check grid dependency, the cases with 24 by 13, 
36 by 20, 48 by 27, 60 by 34, and 72 by 41 nodes are run for jet Reynolds 
number ( Re j ) equal to 600 and Prandtl number equal to 1.25. Figures 2 and 3 
show that there are only minor changes in ( 3 u*/ 3 r*) distribution along the 
wall and (3T*/3x*) distribution across the interface between grid distributions 
60 by 34 and 72 by 41. It is then decided to use 72 by 41 grid nodes for all 
the calculations in the present study. Calculations are performed on a CRAY- 
XMP computer located at NASA Lewis Research Center. The convergent solutions 
are considered to be reached when the absolute value of (eq. ( 2 ) - eq. ( 1 ))/ 
(eq. (1)) is less than 0.003 and the maximum of absolute residual sums for each 
dependent variables is less than 10 ~ 6 . 


RESULTS AND ANALYSIS 

A potential core length, xp, is defined as the distance between the jet 
exit and the point where the jet centerline velocity begins to decrease. As 
shown in figure 4, the potential core length is linearly increasing with jet 
Reynolds number and can be described by the following equation 

_ c Re 
d " C 1 Re j 

where C] is 0.0067 for the tankage considered in the present study. The con- 
clusion of a linear relation between potential core length and jet Reynolds 
number is in agreement with the analysis of free jet in reference 12. The 
centerline velocity, u*^, decays slowly approaching the liquid-vapor interface. 

The velocity drops to zero abruptly at the interface with greater decreasing 
rate for a higher jet Reynolds number. The upstream effect of the existence of 
the interface is less for higher jet Reynolds numbers. It is noted that both 
the potential core length and the decaying rate of centerline velocity are less 
than those in reference 12 because of the difference between two physical flow 
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systems considered resulting in different values of constant C] . As the jet 
approaches the interface, an impingement region is created right underneath the 
interface. This region plays the most important role in the bulk-mixing proc- 
ess since the flow in this region removes the heat from the interface into the 
bulk liquid and induces the vapor condensation. The radial distribution of the 
free surface velocity, v*, is shown in figure 5. The v* increases from zero 

at the centerline to reach a peak shortly and then decreases approaching the 

side wall. For a higher jet Reynolds number, this peak velocity is larger and 

is located closer to the centerline. As expected, the average dimensionless 

surface velocity v* is increasing with Rei and is approximately propor- 

tional to Rel^ 2 : 


where C 2 is 0.0134 in the present study. Consequently, the average interface 
radial velocity is approximately increasing with u?/ 2 . 

The existence of liquid-vapor interface and tank side walls and the with- 
drawing of the liquid from the outer part of the tank bottom enhance the flow 
circulation in the tank. Figure 6 shows the velocity gradient distribution at 
the side wall owing to the liquid flowing toward the tank bottom. Because the 
backward axial velocity is very small the wall shear is negligible at most of 
the tank side walls except the region near the interface. The value of the 
maximum wall shear is higher and is located closer to the interface for jets 
with higher Reynolds number. The distance of the location of maximum wall 
shear from the interface is observed to be approximately the same as the thick- 
ness of the impingement region. From figure 6, the radial boundary layer 
thickness appears to reach a value of about 0.025 D (i.e., about a value of 
0.5 d) as the jet Reynolds number is increased. This finding is in agreement 
with that for a turbulent jet in reference 2. 

For Prandtl number Pr = 1.25, the axial temperature distribution at the 
centerline, R = 0, as a function of jet Reynolds number is shown in figure 7. 
The centerline temperature increases slowly approaching the interface with 
higher increasing rate for a smaller jet Reynolds number. Consequently, in 
the region near the interface, higher jet Reynolds numbers yield much higher 
temperature gradients to approach the interface temperature. The radial dis- 
tribution of temperature gradient at the interface is given in figure 8. As 
observed, a high jet Reynolds number flow induces more interface condensation. 
Most of the mass condensation occurs in the central part of the interface for 
the ratio of jet diameter to tank diameter equal to 0.05 considered in the 
present study. It is expected that if the jet nozzle to tank diameter ratio 
is small enough, the difference of interface configuration between low-g and 
normal-g will not have significant effect on the interface condensation. 
Therefore, as mentioned in the introduction section, the fluid mixing processes 
of low-g and normal-g can be equivalent if the buoyancy force is neglected. 
Calculations were also performed for Prandtl numbers varied from 0.85 to 2.65 
with various jet Reynolds numbers. The effect of Prandtl number on the distri- 
butions of centerline temperature and interface temperature gradient was found 
to be similar to that of jet Reynolds number as shown in figures 7 and 8. 

Higher Prandtl number will yield larger condensation heat flux over the inter- 
face and thinner thickness of thermal layer underneath the interface. It is 
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noted that for higher jet Reynolds numbers and Prandtl numbers, the maximum of 
temperature gradient at the interface is slightly shifted away from the center- 
line due to the effect of significant surface motion. The radial distribution 
of temperature at x* = 0.97 is shown in figure 9. The higher jet Reynolds 
number yields lower flow temperature in the central jet region, approaching 
bulk liquid temperature with higher increasing rate at the edge of jet region, 
and results in a thinner thermal layer underneath the interface. 


The dimensionless condensation-induced velocity at the surface, u*, can 
be calculated by 

u. 


u* = -r- = 


(9T*/9x*> s Ja. 
Pr Re.B 

J 


where the jet Jacob number is defined as Cp(T s ~ ^j^fg' ^he max ^ mum °f u £ 

is located at the centerline. For liquid hydrogen at a saturation pressure of 
30 psia with 2 °R of jet subcooling, the maximum values of u* for Re. = 600 

and 150 are equal to 0.0011 and 0.0009, respectively. From figure 5, it is 
observed that u* is negligibly small compared with surface velocity v* 

except for the very small stagnation region. The temperature difference 
between the bulk liquid and the interface is much smaller than that between the 
jet and the interface. For (T s - Tj) equal to 2 °R, (T s - T ou ^-) is approxi- 
mately equal to 0.06 °R and the corresponding bulk Jacob number, Ja^, is equal 
to 0.001. It can then be justified that the assumption of neglecting of the 
condensation-induced velocity, u c , and bulk Jacob number, Ja^, in the present 
study is acceptable. 


The correlation for the average condensation Nusselt number is assumed to 
be of the form: 


Nu = Re? Pr n (10) 

c 3 j 

Based on the obtained numerical solution, the values of C 3 , m, and n are 
approximately equal to 0.049, 1.0, and 0.97, respectively. As shown in fig- 
ures 10 and 11, the average condensation Nusselt number is linearly propor- 
tional to jet Reynolds number and nearly a linear function of Prandtl number. 

It is noted that the constant n is actually a weak function of Prandtl 
number itself and is decreasing from 1.0 for Pr = 0.85 to 0.95 for Pr = 2.65. 
The value of n is expected to slightly decrease as the Prandtl number keeps 
increasing. Equation (10) also implies that the condensation mass flux is lin- 
early increasing with jet Reynolds number which is in agreement with the data 
for turbulent jet in reference 3. From equations (9) and (10), the correlation 
for average condensation Stanton number can be expressed as 

St^ = C 4 Pr "- 1 = C 4 Pr-0-03 ( 11 ) 

where C 4 is equal to 0.00245. The average condensation Stanton number is a 
weak function of Prandtl number only and is independent of jet Reynolds number. 
In the present study, only the jet Reynolds number and the Prandtl number are 
investigated and the other relevant tankage system parameters such as 
B(= D/d), x s /D, and A ou ^/Aj are specified. Therefore, it is expected that 
the correlation constants, C-| , C 2 , C 3 and C 4 , obtained from the above analysis 
need to be modified if a different tankage system is considered. 
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From equations (2), (5), and (6), the average condensation Stanton number 
for Jat, = 0 can be expressed as 


St 


c 



( 12 ) 


If the outflow temperature T ou t 


approaches T s> equation (12) will become 


St 


c 



(13) 


It implies that St c is nearly a constant and is simply determined by system 
geometry provided that the flow conditions can yield T ou t close to T s . The 
physical phenomenon behind equation (13) is that the cold jet simply enters, 
sweeps by the interface, and receives as much heat as it can absorb, reaching 
saturation temperature and exiting. Equation (13) gives an upper bound for 
St c based on an energy balance with the outflow liquid at its maximum possi- 
ble temperature. In the present computation, d/D = 1/20, and equation (13) 
gives 


St c = 0.0025 (14) 

which is in excellent agreement with the numerical correlation shown in equa- 
tion (11). It can be concluded that for laminar operating conditions with a 
jet Reynolds number up to 600, the injected liquid reaches essentially inter- 
face saturation temperature before exiting (at least for a tankage system of 
d/D = 1/20), and the total condensation rate is controlled by the jet volume 
flow rate rather than by the condensation heat-transfer coefficient at the 
i nterface . 


CONCLUSIONS 

The nondimens ional form of steady state continuity, momentum, and energy 
equations has been solved by a finite-difference method. Numerical solutions 
were obtained for a laminar jet-induced fluid mixing in a cryogenic tank system 
shown in figure 1. The temperature variation in the whole flow field was small 
so that the thermodynamics and transport properties were assumed constant. The 
effect of buoyancy force was neglected. The interface was assumed to be flat 
(wave free) and shear free. The condensation-induced velocity at the interface 
and the associated bulk and jet Jacob numbers were assumed zero. Calculations 
were performed for jet Reynolds numbers ranging from 150 to 600 and Prandtl 
numbers ranging from 0.85 to 2.65. Based on the tankage system considered in 
the present study, it is concluded that: 

(1) The ratio of average radial interface velocity to jet velocity is 
approximately proportional to and the average interface radial velocity 

is approximately proportional to The thickness of the impingement 

J 

region layer underneath the interface is decreasing with increasing jet 
Reynolds number and appears to approach a value of one-half of the jet diameter 
as jet Reynolds number is further increased. 
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(2) If the jet nozzle to tank diameter ratio is small enough, most of the 
vapor condensation will occur at the central part of the interface and the 
effect of the difference of interface configuration between low-g and normal-g 
on the interface condensation will be insignificant. 

(3) The condensation Nusselt number and then the condensation mass flux 
are linearly increasing with jet Reynolds number. The condensation Stanton 
number is independent of jet Reynolds number and is proportional to Pr 1_n 
with n approximately equal to 0.97 for Prandtl numbers ranging from 0.85 to 
2.65. The value of the exponent n is slightly decreasing with increasing 
Prandtl number. The thermal -1 ayer thickness underneath the interface is thin- 
ner for higher values of jet Reynolds number and Prandtl number. 

(4) For laminar operating conditions with a jet Reynolds number up to 600, 
the injected liquid reaches essentially interface saturation temperature before 
exiting (at least for a tankage system with x s /D = 1 and d/D = 1/20), and the 
interface condensation rate is jet volume flow limited. 
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FIGURE 2. - DISTRIBUTION OF DIMENSIONLESS VELOCITY GRADIENT 
AT THE TANK SIDE WALL. 
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FIGURE 3. - RADIAL DISTRIBUTION OF DIMENSIONLESS TEMPERA 
TURE GRADIENT AT THE INTERFACE. 
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FIGURE 4. - DISTRIBUTION OF DIMENSIONLESS AXIAL VELOCITY 
ALONG THE CENTERLINE. 



FIGURE 5. - DISTRIBUTION OF DIMENSIONLESS RADIAL VELOCITY 
AT THE INTERFACE. 
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FIGURE 6. - DISTRIBUTION OF DIMENSIONLESS AXIAL VELOCITY 
GRADIENT ALONG THE TANK SIDE WALL. 
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FIGURE 7.- DISTRIBUTION OF DIMENSIONLESS TEMPERATURE ALONG 
THE CENTERLINE FOR Pr = 1.25. 
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DIMENSIONLESS TEMPERATURE 



FIGURE 8. - RADIAL DISTRIBUTION OF DIMENSIONLESS TEMPERA- 
TURE GRADIENT AT THE INTERFACE FOR Pr = 1.25. 



FIGURE 9. - RADIAL DISTRIBUTION OF TEMPERATURE AT x* = 0.97 
AS A FUNCTION OF REYNOLDS NUMBER FOR Pr = 1.25. 
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FIGURE 10. - AVERAGE CONDENSATION NUSSELT NUMBER AS A FUNC- 
TION OF JET REYNOLDS NUMBER. 



FIGURE 11. - AVERAGE CONDENSATION NUSSELT NUMBER AS FUNC- 
TION OF PRANDTL NUMBER. 
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